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Constraints on the Parameters of the Unseen Pulsar in the 
PWN GO.9+0.1 from Radio, X-Ray, and VHE Gamma-Ray 
Observations 
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Abstract. Radio, X-ray, and H.E.S.S. gamma-ray obser- 
vations of the Galactic Center (GC) composite super- 
nova remnant SNR GO.9+0.1 are used to constrain a 
time-dependent injection model of the downstream elec- 
tron spectrum responsible for the total multiwavelcngth 
spectrum. The effect of spindown power evolution as well 
as nebular field evolution is employed to reproduce the 
present-day multiwavelength spectrum. Assuming a neb- 
ular magnetic field decay model of typical H.E.S.S. -type 
pulsar wind nebulae (PWN), ending with a present-day 
field strength of 6/iG, we obtain an initial spindown power 
of ~ 10 38 ergs/s if we assume a birth period and age of 
43 ms and 6, 500 yr respectively to reproduce the proper- 
ties of the SNR shell. This gives a present-day spindown 
power of ~ 10 37 crgs/s, which agrees well with the present- 
day spindown power derived from X-ray observations. 



\ 1. Introduction 



The H.E.S.S. Collaboration recently reported the 
detection of the composite supernova remnant 
SNR GO.9+0.1 in very high energy (VHE) gamma- 
rays [Aharonian et al. 2005] with a significance of « 13c, 
after 50 hours of observations of the Galactic Center 
(GC) region during March-September 2004 with the 
full telescope system (see also |Aharonian et al. 2006] ). 
The unresolved gamma-ray radiation, which appears 
to be associated with the plerionic remnant core, 
corresponds to a photon flux above 200 GeV of 
(5.7 ± 0.7 stat ± 1.2 sys ) x 10~ 12 /cm 2 /s and luminos- 
ity of ~ 2 x 10 34 ergs/s, making this one of the faintest 
VHE gamma-ray sources. For a radially symmetric 
Gaussian emission region, a 95% confidence limit of 
1.3' ~ 3.2 pc on the source radius is obtained, and 
similarly for emission from a uniform thin shell, with 
extension < 2.2' ~ 5.4 pc. This is the first time that 
SNR GO.9+0.1 has been detected at gamma-ray en- 
ergies pQ, subsequent to the earlier unconstraining 
upper limit given by HEGRA [Aharonian et al. 2002] . 
and H.E.S.S. observations in 2003 (resulting in a 4a 
detection) [Aharonian et al. 20"04] . MAGIC also detected 



a small excess from GO.9+0.1 [Albert et al. 2006] . but it 
is however not statistically significant yet due to limited 
observation time. 

Helfand and Becker [5] observed GO.9+0.1 in 1984 for 
45 minute integrations at 20 cm and 6 cm with the VLA 
in their search for further examples of composite objects, 
despite its omission from Green's first catalogue [7] . This 
led to the discovery of the composite nature of this bright, 
extended source near the GC in the radio band, with the 
flat radio core, corresponding to the pulsar wind nebula 
(PWN) , and steeper shell components clearly distinguish- 
able. (Shortly afterwards, Reich, Sofue & Fuerst [12] con- 
firmed G0.9+0.1's non-thermal nature and its identifica- 
tion as an SNR by detailed comparison of radio continuum 
and infrared emission from the GC). Using various previ- 
ous observations, Helfand and Becker assigned spectral 
indices of a c ~ 0.1 (typical of Crab- like remnants) and 
a s k 0.6 to the core and shell components, and also con- 
structed a broad band spectrum (radio - X-ray) of this 
source, as well as an upper limit of 3 Jy on the 25fim 
infrared flux. (Recent VLA measurements at 90 cm [9] re- 
vealed similar indices of 0.12 and 0.77 for the core and 
shell components, confirming this SNR's composite clas- 
sification [55]; see [3] for a recent review). 

Helfand and Becker [5] furthermore indicated a core- 
to-shell diameter ratio of ~ 0.25. Assuming a distance of 
10 kpc, they obtained radio luminosities of £_r. c = 7.0 x 
10 34 crgs/s and L R<S = 2.8 x 10 34 ergs/s for the core and 
shell components, and an X-ray luminosity (0.5 — 4.5 kcV) 
of ~ (1 — 5) x 10 34 ergs/s for this remnant, setting the X- 
ray-to-radio flux ratio for the core to < 0.6. Other cited lu- 
minosity ratios are Lr^ c /Lr iS ~ 2.5 and Lx, c / Lx,s ~ 0.5. 
These authors also argued for an initial spin period of 
a few ms for the putative pulsar, but further VLA ob- 
servations at 6 cm and 20 cm did not reveal any point 
source above 6a (0.45 mJy) [J]. Observations during a con- 
tinuum GC survey at 843 MHz by MOST [S] confirmed 
a core and shell component, with respective flux densities 
in close agreement with the spectra given by [5] . Owing to 
G0.9+0.1's structure being similar to that of SNR 0540- 
693 (in the Large Magellanic Cloud) which contains a pul- 
sar, GO.9+0.1 was identified as a likely candidate for being 
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a pulsar host. High dispersion measures and interstellar 
scattering may however inhibit detection of pulsed emis- 
sion from GO.9+0.1 0. 

While performing a survey of the GC, BeppoSAX 
serendipitously detected SNR GO.9+0.1 Q3]. This detec- 
tion represented the first firm evidence of X-ray emission 
from GO.9+0.1, in addition to an earlier marginal detec- 
tion by the Einstein Observatory [H] and a nondetection 
by ROSAT [TDHH]. Follow-up observations of the GC and 
on-axis observations of GO.9+0.1 [H [23 [23l [24] confirmed 
the non-thermal nature of X-rays from its radio core, and 
thus the composite nature of this remnant, and also found 
marginal indications of a diffuse character. Only an upper 
limit for X-ray emission from the SNR shell could be given, 
probably due to line-of-sight absorption [16] . Although no 
flux variations or pulsations were found, Sidoli et al. [24j 
argued for the existence of a young central neutron star 
based on the broad band (radio & X-ray) properties of the 
SNR. 

The ASCA GC survey also detected significant X-ray 
emission from GO.9+0.1 in the 3 — 10 kcV band, but no 
significant emission in the softer energy band [18] . The 
emitting region was compact, in accordance with the radio 
core, but unresolved, and no line-like feature was observed. 
Similar updated results for the 0.7 — 10 keV band, as well 
as an upper limit of 1.5' for the apparent X-ray size, are 
given in [15] , 

A higher resolution 35 ksec observation of this source 
in October 2000 with Chandra X-Ray Observatory [5] 
revealed a clear axial symmetry for the PWN, which 
matched the morphology of the radio nebula, as well 
as a faint X-ray point source along the symmetry axis 
(CXOU J174722. 8-280915) which is considered to be the 
best candidate for emission from a central pulsar (al- 
though no pulsations were detected). A bright clliplical 
5" x 8" clump, which may be an intermediate-latitude fea- 
ture in the pulsar wind, was also observed. An equatorial 
torus and axial jet morphology provides a natural expla- 
nation for the data and lends support to the idea that this 
type of morphology might be ubiquitous within the pulsar 
population. 

SNR GO.9+0.1 was also easily identified in the hard en- 
ergy band when a wide-angle XMM-Newton survey of the 
CG Region was conducted during 2000-2001 (11 pointings 
of 10 — 25 ksec exposures) [El [17] . On-axis observations 
during September 2000 (~ 29 ksec) [11] revealed a large 
scale X-ray morphology in good correspondence with the 
20 cm VLA radio contours [8] , with the X-ray core match- 
ing the eastern radio peak, and the western side of the 
X-ray arc-like feature corresponding to the western radio 
peak. Diffuse emission from the radio shell was also de- 
tected in X-rays for the first time (of cither thermal or 
non-thermal origin), in good agreement with the earlier 
upper limit of Sidoli et al. [24] , Both an absorbed power- 
law and a thermal bremsstrahlung model provide good fits 
for the data in the region of the PWN. Evidence for spec- 



tral steepening with increasing radius supports a scenario 
where high-energy electrons undergo synchrotron losses as 
they diffuse through the nebula. The spectral steepening 
is stronger than in other PWN (e.g. 3C58 & G21. 5-0.9) 
and the core is harder. Small-scale structure resembled 
that detected by Chandra, having less detail but double 
the extent of the emission. Spectra of the various struc- 
tures were derived for the first time, indicating that the 
eastern part of the arc-like feature exhibited harder X-ray 
emission than the western part. This may be explained in 
terms of relativistic beaming or Doppler boosting effects. 
The jet-like feature also had distinct spectral indices for its 
northern and southern part. Furthermore, the central pul- 
sar candidate CXOU J174722. 8-280915 was visible above 
6 keV, and magnetospheric emission seems to be preferred 
to a black-body model, implying a luminosity of ~ 10 33 d^ 
ergs/s for this source (d\o = d/10 kpc). The deepest X- 
ray observation ever performed on this SNR [22] confirms 
the spectral softening at larger radial distances from the 
PWN peak as well as X-ray emission from a region spa- 
tially coincident with the radio shell. 

During 2003-2004, INTEGRAL observed the GC in 
the 20 — 100 keV energy range for an effective exposure 
time of 4.7 x 10 6 sec [Belanger et al. 2006| . Emission from 
a region with centroid label IGR J17475-2822, a source 
associated with Sgr B2, was clearly observed [13]. The 
extension of this source toward the north was tentatively 
associated with GO.9+0.1 and interpreted as a detection 
of VHE synchrotron radiation. 

In this paper we use radio, X-ray, and gamma-ray data 
to constrain both a time-dependent injection model of the 
downstream electron spectrum (responsible for the total 
multiwavclcngth spectrum), and the evolutionary history 
of the spindown power. 

2. The Model 

We assume that the central pulsar's spindown power has 
the following functional form: 



L = Iuilo = —Kuj 



71+1 



(1) 



with / the moment of inertia, lo the angular frequency, 
oj its time-derivative, and n = ujuj/lo 2 the braking in- 
dex. We next assume that n is a constant (see, however, 
[Ch en fc Li 2006] ) and that the neutron star's crustal 
magnetic field does not decay. This leads to the condition 



PP 



n — 2 



n-2 



(2) 



with P = 2it/u! the pulsar period, P its time-derivative, 
and the subscript '0' indicating quantities at pulsar birth. 
Upon integration of u> from eq. ([1]) , and using eq. ([2]) , one 
finds the following general expression for the evolution of 
L(t) (see also fi"4]): 



Lit) = L 



(n - l)P 2 L t 



An 2 1 



(3) 
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Model Parameter 


Symbol 


Value /Range 




4C- 


Braking index 


n 


3 




: 


B-field parameter 


a 


0.5 


J? 




Present-day B-field 


B(T) 


6 fj,G 




Conversion efficiency 


e 


0.1 




?ni 


Age 


T 


6,500 yr 


_ f 




Characteristic time scale 


To 


500 yr 




38 Ji 1 


Distance 


d 


8.5 kpc 


□ 




Magnetization parameter 


a 


0.2 


£ 

3 


- 


Moment of inertia 


I 


10 45 g.cm 2 


^1 
C 




Q break energy 


E h 


10~ 8 - 10 ergs 




37-5.- 


Q index 1 


ai 


-1.0 


| 




Q index 2 


a.2 


-2.2 


a. 
us 




Initial spindown power 
Birth period 


Lo 
Po 


10 36 - 10 40 ergs/s 
0.043 s 


s 


36 .5 


Table 1. Values of parameters used for the model (see 




text for details). 











In what follows, we assume Po = 0.043 s, as inferred by 
Van der Swaluw and Wu [25] using the ratio of the PWN 
and shell radii. We model the effect of a time-changing 
nebular field by 



B(t) 



Bo 



1 + (t/T ) a 



(4) 



For the electron injection spectrum Q(E,t) (number of 
electrons per second per energy at the shock radius r$), 
we assume a broken power law with indices a\ and a2, 
and break energy E^: 



Q(E,t) 



E<E h 
E>E h 



We normalise Q by requiring that 

J Q(E,t)EdE = eL(t), 



(5) 



with e ~ 0.1 a conversion efficiency of spindown power 
into particle luminosity. This leads to 



Qo(t) 



{at + 2) (q- 2 + 2) 
(a 2 - ai) 



(6) 



for a2 < —2. 

We next impose two boundary conditions. The first is 
that rL<0.5rs, with the electron's Larmor radius [2]. 
This ensures particle confinement within the PWN, and 
leads to a condition on the particle energy: 



E(t)< 



aL(t) 



2 V (l + a)c 



(7) 



with a = Leb / Le the ratio of the electromagnetic energy 
flux to the particle energy flux. The second energy con- 
dition is required to ensure that the particles radiating 
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-a -4 -2 

Break Energy (ergs] 

Fig. 1. Plot of log 10 Xtot (with \tot the rms-total of the 
individual test statistics for fits to the radio, X-ray and 
gamma-ray data) vs. log 10 of initial spindown luminosity 
Lq and log 10 of injection spectrum break energy E^, in 
order to find the 'best fit' (indicated by an 'X') for Lq and 
£"b (at a fixed age T = 6,500 yr and birth period Po — 
0.043 s). 



synchrotron emission will survive until time t = T (with 
T the PWN's age): 



E(t)< 



422 



B{t)B{T){T - t) 



ergs. 



(8) 



We lastly calculate the leptonic spectrum dN/dE subject 
to eq. ([7]) and eq. §8§ by integrating Q over time from 
t = to t = T = 6,500 yr (see [TO], where an age of 
~6,800 yr is derived assuming that G09.+0.1 is expanding 
adiabatically). Table [TJ summarizes the values adopted 
for certain model parameters. 



3. Results and Conclusions 

Using this time-dependent model for the nebular aver- 
aged leptonic spectrum, we add the contributions from all 
epochs which survive to the present-day to give the net 
present-day photon spectra as calculated for synchrotron 
and inverse Compton emission. In the latter case the tar- 
get photon fields are the CMBR, as well as galactic tar- 
get photon fields from the 25K dust and starlight photon 
fields, assuming associated energy densities of ~ 1 eV/cm 3 
for each component. 

From Figure [TJ we find an optimum fit of initial spin- 
down luminosity Lq ~ 2.5 x 10 38 crgs/s and a spectral 
break energy ofPb~8xlO~ 3 ergs. Note that we ob- 
tain different optimum values when we fit the radio, X- 
ray and VHE gamma-ray data individually. The resulting 
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All 



Gamma-rays 




Energy (eV) 

Fig. 2. Plots of best fits of the spectral energy distribution 
(SED) of GO. 9+0.1, with the different fits corresponding 
to Table [2] i.e. for individual and collective fits of the 
multiwavelength data (at a fixed age T = 6,500 yr and 
birth period Po = 0.043 s). Also shown are radio data 
points (the diamonds) and radio spectrum taken from [3], 
X-ray spectrum from [22| , and gamma-ray spectrum from 
[Aharonian et al. 2005| . 

fits (and multiwavelength observations) are shown in Fig- 
ure [2] and Table [2] where we also show the correspond- 
ing total energy output Ftot = / L(t)dt, integral X-ray 
flux Fx (2 — 10 keV) in crgs/cm 2 /s, and integral 7-ray 
flux above 0.2 TeV, F 7 (> 0.2), in /cm 2 /s. A distance 
of d = 8.5 kpc is assumed [Aharonian et al. 2005| . These 
values correspond quite well to inferred values given in 
[Aharonian et al. 20051 fTUl 121 HT1 123] 1. 

Eq. ([3]) now implies a value of present-day spindown 
power of L(T) ~ 7.4 x 10 36 crgs/s, which agrees within 
~ 50% with the inferred value of L(T) ~ 1.5 x 10 37 ergs/s 
from BeppoSAX observations [24j . One may then also infer 
a pulsar polar cap (PC) magnetic field strength of Bpc ~ 
3.2 x 10 1! VPoPo ~ 5 x 10 12 G, which is a typical value 
for the canonical pulsar population. 

We are, however, not yet satisfied with the quality of 
the total fit, since it is obvious that the radio data are not 
fitted very well. This is also the case for an alternative 
nebular field which is initially larger, and falls exponen- 
tially to the present-day assumed value of B(T) = 6fiG. 
Future investigations will relax the initial spin period Po 
and age T (keeping them free), and include extra com- 
ponents in the injection spectrum, to see if we can ob- 
tain better quality fits for the spectral energy distribution 
(SED) of GO.9+0.1. 

Acknowledgements. This publication is based upon work sup- 
ported by the South African National Research Foundation 
under Grant number 2053475. 
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Data 


E h 


L 


F 7 (> 0.2) 


Fx 


F to t 


X 


7-rays 


0.5 


38.0 


-11.4 




48.8 


-1.0 


X-rays 


1.0 


37.7 




-11.6 


48.7 


-0.35 


Radio 


-2.1 


39.8 






49.0 


2.2 


All 


-2.1 


38.4 


-11.9 


-12.4 


48.9 


2.4 


Obs. 






-11.2 


-11.2 







Table 2. Table of 'best fit' parameters for individual and 
collective multiwavelength fits, as well as observed val- 
ues [Aharoni an et al. 20051 [TT] . Column units are log 10 of 
ergs, ergs/s, /cm 2 /s, ergs/cm 2 /s, ergs, and log 10 of a di- 
mensionless test statistic. (See text for details). 
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